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ABSTRACT 

A  one-dimensional  coupled  chemistry-radiation  model  (ChemRad)  for  predicting 
combustion  and  heat  transfer  characteristics  in  a  single  and  multi-layered  porous 
media  is  presented.  The  model  allows  detailed  gas  and  surface  chemical  kinetics 
mechanisms  to  be  incorporated,  and  uses  a  two-flux  radiation  approximation  in  the 
energy  equation.  A  stoichiometric  propane/ air  mixture  burning  in  a  thin  wire-mesh 
burner  has  been  modelled  using  a  well-documented  high-temperature  propane 
kinetics  scheme  with  40  reactive  species  and  93  reactions.  A  systematic  analysis  of  the 
solution's  sensitivity  to  variations  in  input  parameters,  particularly  the  optical 
characteristics  of  the  media,  has  been  performed  and  influential  parameters  have  been 
identified.  It  has  been  found  that  the  dominant  factors  that  influence  the  solution,  in 
order  of  significance  are  the  area  density,  voidage  and  reflectivity  of  the  porous 
material,  and  the  ratio  of  forward-to-backward  scattering  of  radiation. 
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A  Parametric  Analysis  of  a  Coupled  Chemistry- 
Radiation  Model  in  Porous  Media 


Executive  Summary 

The  accuracy  and  reliability  of  numerical  models  to  predict  combustion  and  thermal 
radiation  characteristics  of  gaseous  and  solid  fuels,  relies  heavily  on  their  ability  to 
adequately  couple  the  various  transport  mechanisms  in  the  process.  These  processes 
usually  include  chemical  reactions,  turbulent  mixing,  and  radiation  heat  transfer.  The 
work  presented  here  is  aimed  at  developing  a  numerical  model  ( ChemRad)  to  calculate 
the  laminar  combustion  inside  porous  media,  accounting  for  detailed  chemistry  and 
radiative  heat  transfer.  This  study  is  part  of  an  ongoing  research  program  aiming  at 
developing  computational  capabilities  to  model  the  performance  of  infrared  (IR) 
emitting  devices.  More  specifically  to  assess  the  potential  use  of  porous  radiant  burner 
technologies  in  the  design  of  a  new  generation  of  IR  tracking  and  decoy  devices. 

The  numerical  model  presented  here  has  a  unique  capability  for  handling  complex  and 
detailed  chemical  kinetics  accurately  and  efficiently,  hence  eliminating  chemistry- 
related  errors  in  the  solution.  The  greatest  source  of  errors  in  the  model,  however 
remains  the  uncertainty  in  some  input  parameters  that  describe  the  microstructure  and 
optical  properties  of  the  porous  medium.  Accordingly  a  systematic  analysis  of  the 
solution's  sensitivity  (or  lack  of  it)  to  variations  in  nine  input  parameters  has  been 
conducted.  Three  of  the  parameters;  voidage,  tortuosity,  and  area  density  represent 
the  microstructural  properties  of  the  porous  medium,  and  the  remaining  six 
parameters  describe  the  optical  properties  of  the  porous  material.  These  are  the 
emissivity,  reflectivity,  single  scattering  albedo,  extinction  coefficient,  forward-  and 
backward  radiation  scattering  fractions.  It  has  been  found  that  the  most  influential 
parameters  (in  order  of  significance)  are  the  area  density,  voidage,  reflectivity,  and  the 
forward-to-backward  scattering  ratio. 

The  capability  of  the  ChemRad  model  to  capture  the  combustion  characteristics  has 
been  successfully  demonstrated  by  a  favourable  agreement  between  the  model 
predictions  and  the  experimental  results,  which  are  documented  in  a  separate  technical 
report  to  be  published  soon. 
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1.  Introduction 

During  the  last  decade  the  utilisation  of  Porous  Radiant  Burner  (PRB)  technology  has 
emerged  as  a  practical  method  for  developing  combustion  systems  with  low  hazardous 
emissions  gases,  such  as  CO  and  NOx  [1-3].  Another  advantage  of  a  PRB  is  its  capability 
to  maintain  stable  combustion  of  premixed  flame  and  enhanced  radiation  heat  transfer 
over  a  wider  operating  range  than  a  conventional  "blue-flame"  premixed  burner  [4]. 
PRB  technologies  have  been  primarily  driven  by  the  need  to  improve  the  thermal 
efficiency  of  existing  heat  and  power  generating  facilities  and  also  by  the  restrictions  on 
the  emission  levels  of  hazardous  gases.  To  the  best  knowledge  of  the  author  no 
published  research  exist  that  attempts  to  examine  the  potential  use  of  PRB  technology 
for  Defence  applications,  for  example  as  Infrared  (IR)  decoys  or  tracking  devices.  The 
"Hot  Brick"  technology  uses  a  similar  concept  to  the  PRB,  but  it  electrically  heats  a 
solid  block  of  metal  to  a  high  temperature  causing  it  to  emit  in  the  IR  wave  band.  It  is 
anticipated  that  the  present  work  will  deliver  numerical  capability  that  will  assist  in 
determining  the  suitability  of  PRB  technology  for  applications  as  IR  emitting  devices. 

The  objectives  of  this  study  are  firstly  to  develop  a  numerical  model  to  simulate  the 
combustion  process  in  porous  media  accounting  for  coupled  chemistry-radiation 
processes.  Secondly  to  perform  a  sensitivity  analysis  of  the  model's  parameters, 
particularly  those  that  are  difficult  to  accurately  determine  experimentally  or 
analytically.  The  capability  and  robustness  of  the  model  will  be  verified  with 
measurements  taken  from  an  experimental  PRB.  This  however  is  reported  in  a  separate 
document  [5].  This  work  was  conducted  as  part  of  an  AIR  sponsored  task,  AIR  97/250 
"Modelling  IR  Countermeasures  for  RAAF". 


2.  PRB  Concept 

Porous  radiant  burners  are  manufactured  either  from  ceramic  or  stainless  steel  fibers  or 
from  fired  clay.  Approximately  90  to  95  percent  of  the  porous  medium  consist  of  micro 
pores.  The  porous  medium  being  typically  at  around  900  degrees  centigrade  radiates 
energy  and  depending  on  the  emissivity  of  the  material,  heat  flux  intensity  of  more 
than  300kW  per  square  meter  can  be  achieved.  In  principle,  the  combustion  process  in 
a  PRB  is  a  premixed-based  combustion  but  with  a  fundamental  difference  from  a  "blue- 
flame"  burner  in  the  mechanism  in  which  heat  is  transferred.  Unlike  conventional 
premixed  flame,  shown  schematically  in  Fig.l  (top),  in  which  heat  convection  is  the 
dominant  mode  of  heat  transfer,  PRB  transfers  heat  mainly  by  thermal  radiation  as 
illustrated  in  Fig.l  (bottom).  Close-up  photographs  of  blue-  and  porous-burner  flames 
are  shown  in  Fig.  2,  and  clearly  explain  the  origin  of  the  term  'blue  flame'. 
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Figure  1:  Schematic  Illustrations  of 'blue  flame'  (top)  and  porous  radiant  flame  (bottom). 

In  PRBs  some  of  the  radiant  heat  is  also  recirculated  to  the  unbumt  mixture.  Such 
recirculation  of  heat  without  simultaneous  dilution  by  the  burnt  gases  has  the  effect  of 
extending  the  flammability  limits  and  improving  flame  stability  [6-8].  The  heat  transfer 
from  the  reaction  zone  inside  the  porous  medium  to  the  solids  causes  the  solids  to 
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radiate  thus  lowering  the  flame  temperature  and  consequently  reducing  (thermal)  NOx 
emissions. 

The  improved  characteristics  of  PRB  combustion  include  increasing  the  flame  speed, 
decreasing  the  lean  flammability  limit  and  enabling  the  burning  of  low  heat-content 
fuels.  A  number  of  studies  [9,10]  showed  that  CO  is  created  and  "consumed" 
(converted  into  CO2)  within  the  porous  medium.  They  also  showed  those  important 
radicals  such  as  OH,  H  and  O,  peak  within  the  porous  layer  and  almost  all  NO 
formation  occurs  within  the  porous  medium  because  the  gas  temperature  loses  enough 
heat  to  the  solids  to  suppress  the  thermal  NOx  reactions.  To  illustrate  the  effect  of 
thermal  radiation  on  the  combustion  temperature,  a  comparison  of  gas  temperature  in 
a  porous  burner  and  a  conventional  burner-stabilised  blue  flame  is  shown  in  Fig.  3. 
These  calculations  are  carried  out  using  the  ChemRad  model,  which  will  be  described  in 
a  later  section,  and  refer  to  a  stoichiometric  propane/ air  mixture  and  identical  inlet 
conditions  ( e.g .  reactants  composition,  temperature,  flow  rate  etc.)  for  both  burners.  It  is 
clear  from  the  figure  that  the  flame  temperature  in  a  porous  burner  is  lower  than  that 
in  a  'blue-flame'  burner  by  approximately  700K.  This  difference  is  primarily  attributed 
to  the  enhanced  radiative  heat  transfer  of  the  porous  burner. 


Figure  2.  A  close-up  photograph  of  Sydney  University's  "Blue  Flame"  burner  (left)  and 
DSTO's  Porous  Radiant  Burner  (right). 


3 


DSTO-RR-0188 


Figure  3:  A  comparison  between  the  adiabatic  gas  temperature  in  a  porous  burner  and 

that  in  a  burner-stabilised  blue  flame,  for  a  stoichiometric  propane/air  mixture. 


3.  The  Numerical  Model 


The  numerical  model  is  one-dimensional,  laminar,  adiabatic  model.  It  includes  a 
separate  energy  equation  for  the  gas  and  the  solid-phase  accounting  for  conduction, 
convection,  and  radiation  heat  transfer  mechanisms  in  all  phases.  The  model  allows 
detailed  chemistry  of  gaseous  species  and  surface  reactions  to  be  incorporated. 


3.1  Model  Description 

The  coupled  chemistry-radiation  model  ChemRad,  is  capable  of  handling  a  single  or 
multi-layered  burner  of  different  porous  media  and  applies  for  planar,  cylindrical  or 
spherical  geometry.  The  ChemRad  model  is  based  on  the  work  of  Shardlow[ll],  and 
utilises  the  core  routines  of  the  chemical  kinetics  package,  CHEMKIN  [12-15].  In  the 
present  study  the  porous  medium  has  been  represented  as  a  thin  three-layered  mesh  of 
the  same  material.  The  geometry  and  material  used  in  the  model  have  been  selected  to 
match  those  of  an  experimental  PRB,  which  has  been  investigated  in  a  separate  study 
[5],  A  schematic  description  of  the  model  is  shown  in  Fig.  4. 
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Figure  4:  Schematic  description  of  the  computational  domain  for  the  porous  burner  model. 


Users  are  required  to  specify  the  reactants  compositions,  inlet  temperature,  burner 
(computational)  dimensions,  and  the  physical  and  optical  properties  of  the  porous 
material. 


3.2  Governing  Equations 

The  governing  equations  for  mass,  species,  and  energy  conservation  are  described  next. 
In  these  equations,  p  is  the  gas  density,  u  is  the  gas  velocity,  <J>  is  the  porosity  of  the 
solid  (<(>=1  in  the  gas  region),  Yk,  Vk ,  TV*  and  ox  are  the  mass  fraction,  diffusive  velocity, 
molecular  weight,  and  the  production  rate  of  the  kth  species,  respectively.  It  should  be 
mentioned  here  that  the  governing  equations  presented  below  do  not  include  the 
surface  reactions  terms,  which  have  not  been  accounted  for  in  this  study.  A 
comprehensive  description  of  the  governing  equations,  including  the  terms  for  surface 
chemistry  can  be  found  in  Ref.  [11]. 

3.2.1  mass  continuity  equation 

m"=  pu(j> ,  (1) 

where  m" ,  is  the  mass  flow  rate  per  unit  area. 

3.2.2  gas  species  conservation  equation 

pu^^f.+4-(purlvt )  =  ®,  wj...k  =  m  (2) 

ax  ax 

where  kk,  is  the  total  number  of  gaseous  species. 
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3.2.3  gas-phase  energy  equation 


pu(f>c 


PS 


dl±. 

dx 


d_ 

dx 


f 


dx 


dT  \  “  dT  « 

8  +  J ]p<t>YtVtcrs-f-  +  h,(Ts-T,)  =  -^mt  htWJ,(3) 
J  k= 1  “-x  k=\ 


where  cpg  is  the  specific  heat  of  the  kth  gaseous  species,  Tg  ,  Ts  are  the  gas  and  solid 
temperatures,  respectively,  hv  is  the  volumetric  heat  transfer  coefficient,  hk  is  the 
enthalpy  of  the  kth  species  and  ks,e  is  the  effective  thermal  conductivity  of  the  gas 

(sw^/Cg). 

The  production  rate  term  (&*)  is  calculated  using  the  kinetic  mechanism  that  describes 
the  chemistry  of  the  fuel  .  A  general  form  for  representing  the  chemistry  of  Wc* 
chemical  species  is: 

kk  kk 

IX,.  A,  ^  IX  A  * ---O'  =  h-J)  /  (4) 

k= 1  k= 1 


where  are  the  stoichiometric  coefficients,  I  is  the  total  number  of  reactions,  and  A k  is 
the  chemical  symbol  of  the  kth  species.  The  superscripts  '  and  "  refer  to  the  forward 
and  backward  stoichiometric  coefficients,  respectively.  Accordingly,  the  term  for  the 
production  rate  for  the  set  of  chemical  reactions  described  by  Eq.  (4),  is  given  as: 


C Ok 


kk 

kf£[[xtr 


k= 1 


kk 

-*«n  [**]■'■ 


k= 1 


■■■{k  =1  ,...,kk), 


(5) 


where  [Xk]  is  the  molar  concentration  of  the  kth  species,  and  kfi ,  hi  are  the  forward  and 
backward  rate  constants  of  the  ith  species,  respectively.  The  forward  rate  constant,  k/  is 
calculated  using  the  Arrhenius  temperature  dependence  formula: 


kn-  yf,.rA  exp! 


-Ej 

RT 


(6) 


In  Eq.  (6)  the  pre-exponential  factor  A,,  the  temperature  exponent  /?, ,  and  the  activation 
energy  E, ,  are  specified  by  the  kinetic  mechanism.  The  backward  rate  constant  hi,  is 
calculated  indirectly  from  the  equilibrium  rate  constant  hi : 

K  =  if-’  <7> 

Ki 

where  the  equilibrium  constants  h  is  given  by: 


MKi)  =  Z (°'ki  -  °ki)AG't‘  ~ (v'ki  -  v"ki )  ln(^)  • 

i=i  Ai 


(8) 
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The  Gibbs  free  energy  AGh*  (normalised  by  RT,  R  being  the  universal  gas  constant)  is 
calculated  from  the  changes  in  the  system's  normalised  enthalpy  AHG  ,and  entropy 
ASu-  These  terms  are  correlated  in  the  relationship: 

AG’k,  ~  AH’ id  -  TAS*ki .  (9) 

The  normalised  enthalpy  and  entropy  terms  are  calculated  using  the  polynomial 
formulation: 

a +^-  an) 


A5«=ZfrV+aI.> n(r.)+"™ 

ti  (»- 1) 


where  the  polynomial  coefficients,  an,  are  obtainable  from  a  number  of  databases  [16- 
18], 

3.2.4  solid-phase  energy  equation 


4-\K.^r\+K(Ts-T.)=^r 

ax  ax  ax 


The  terms  in  Eq.  (12)  represent  conduction,  convection  and  radiation  heat  transfer, 
respectively.  An  empirical  correlation  has  been  used  to  determine  the  convective  heat 
transfer  coefficient  hv  ,  for  the  porous  medium.  In  the  current  model  of  a  metal  fibre 
matrix,  the  Nusselt  number  has  been  estimated  using  the  following  correlation[19]: 


Nud  =0.10  Re/64, 


where  the  Nusselt  and  Reynolds  numbers  are  calculated  based  on  the  average  fibre 
diameter  and  room  temperature.  The  effective  thermal  conductivity  of  the  solid  ke,s  is 
evaluated  by  the  following  correlation  [11]: 


—  =  1  + - - - , 


where  the  parameter  m,  is  given  by: 


1  1 

m=  1.20-29—  (0.81  -<t>)2  +  1.09  +  2.5 — . 

Lf  Lf 
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In  Eq.  (15),  d  and  L/,  are  the  diameter  and  the  length  of  the  fibre. 

3.3  Radiation  Model 

A  two-flux  radiation  model,  also  known  as  the  Schuster-Schwarzchild  approximation 
has  been  used  to  calculate  the  net  radiative  flux  qr,  in  Eq.  (12).  The  model  is  valid  for 
radiative  transfer  in  absorbing,  emitting  and  scattering  media.  The  main  limiting 
assumption  of  the  model  is  that  the  medium  properties  are  such  as  to  yield 
homogenous  absorbing,  emitting  and  scattering  of  electromagnetic  radiation  in  a 
forward  scattering  hemisphere  and  in  a  backward  scattering  hemisphere,  separately. 

The  net  radiative  flux  qr/  is  calculated  as  the  difference  between  the  forward  radiation 
flux  q+,  and  the  backward  radiation  flux  q~,  (qr  =q+  -  q~).  These  radiative  transfer  fluxes 
are  obtained  by  solving  the  following  equations: 

”7  =  _2cr^+  +  2(7 eao  (fq  +bq~)  +  2cre  (1  -  a>0  )oTf4  (16) 

and, 

=  ~2aeq~  +  2 cyeco0(fq-  +  bq+)  +  2cre(l -o)0)aT* .  (17) 

Here  coo  is  the  single  scattering  albedo,  which  represents  the  fraction  of  attenuated 
energy  that  is  the  result  of  scattering.  For  a  non-scattering  medium  coo=0,  while  for  pure 
scattering  coo  =1.  ae  ,  is  the  extinction  coefficient,  which  represents  the  fraction  of 
attenuated  energy  that  is  the  result  of  absorption  and  scattering,  and  f  &  b  are  the 
forward-  and  backward-radiation  scattering  fractions,  respectively,  a,  is  the  Stefan- 
Boltzmann  constant. 


3.4  Boundary  Conditions 

The  boundary  conditions  required  for  solving  the  governing  equations  consist  of  user 
specified  inlet  conditions  and  a  set  of  conditions  dictated  from  energy  balance 
considerations.  At  the  inlet  (x=xin)  the  gas  temperature  and  the  mass  fraction  of  the 
reactants  are  known  quantities:  Tg=Tg,in/  Yk=Yk,in .  The  boundary  condition  at  the 
upstream  surface  of  the  porous  layer  x=xo  ,is  established  from  an  energy  balance 
equation: 


(1  -  m,  ^  L„,  =  (1  -  (Tt  -  T, )  - [q- ■  (x, )  - q' (x0 )] , 


(18) 
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where  ho,  is  the  convective  heat  transfer  coefficient  per  unit  area  and  ks,  is  the  thermal 
conductivity  of  the  solid.  The  radiative  heat  transfer  components  are  calculated  using: 

q+  O0)  =  scT*  +  XS  (xo)  (19) 

where  T.  is  the  upstream  temperature  of  the  environment  to  which  the  upstream 
boundary  of  the  porous  medium  is  exposed,  i.e.  burner  housing,  yo  ,is  the  reflectivity  of 
the  porous  medium,  and  s,  is  the  average  emissivity  of  the  burner  surface,  q-(xo)  in  Eq. 
(17)  is  the  backward  radiation  fraction  at  xo. 

The  boundary  condition  at  the  downstream  surface  of  the  burner  (x=xi)  : 

q~(xL)  =  eaT*  +  x0q+  (xL ) ,  (20) 

T+  being  the  downstream  temperature  of  the  environment  to  which  the  downstream 
boundary  of  the  porous  medium  is  exposed  i.e.  the  ambient  temperature.  At  the  outlet 
x=xout,  the  boundary  conditions  are  dTg/dx=0  and  dYi/dx=0.  These  conditions  imply  that 
the  chemical  reactions  are  completed  and  the  process  is  adiabatic. 


3.5  Numerical  Solver 

The  method  for  solving  Eqs.  1-20,  is  to  discretise  the  conservation  equations 
establishing  a  set  of  simultaneous  algebraic  equations  that  are  applied  on  discrete  mesh 
points  within  a  predetermined  computational  domain.  The  core  solver  is  based  on 
Sandia's  PREMIX  module[14],  and  makes  use  of  the  well-known  CHEMKIN 
package  [12-15].  It  uses  a  combined  time-dependent  and  steady  state  method.  That  is,  a 
modified  damped  Newton  method  is  initially  used  to  solve  the  non-linear  algebraic 
equations  by  an  iterative  process.  If  the  Newton  method  fails  to  achieve  the  required 
level  of  convergence,  then  a  time-stepping  algorithm  is  used.  The  idea  is  to  determine 
the  steady-state  solution  by  solving  transient  equations,  namely  for  each  of  the 
conservation  equations  (mass,  species,  etc..)  an  additional  time-dependent  component 
is  added.  For  example,  an  extra  term  dYi/dt  is  added  in  the  left-hand  side  of  Eq.  2.  The 
combination  of  Newton  and  time-stepping  methods  utilises  the  advantage  of  the  rapid 
convergence  of  Newton's  method,  and  the  robustness  of  the  time-stepping  algorithm 
for  handling  the  stiffness  problem. 

An  adaptive  meshing  placement  is  a  powerful  tool  to  alleviate  convergence  difficulty 
and  reduce  the  user  intervention  in  influencing  the  convergence  of  the  solution.  The 
iterative  process  commences  with  a  coarse  mesh,  which  is  easy  to  converge  and 
provides  an  improved  guess  over  that  specified  by  the  user.  The  refinement  of  mesh  is 
then  done  automatically  allowing  more  mesh  points  to  be  inserted  in  regions  with  high 
gradients  and/  or  curvatures  in  the  concentration  profiles.  The  gradient  (GRAD)  and 
curvature  (CURV)  parameters  have  a  strong  effect  on  the  convergence  rate  (but  not  on 
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the  solution's  accuracy).  The  parameter  GRAD  determines  the  number  of  additional 
mesh  point  to  be  inserted  in  regions  with  high  gradients  of  concentrations,  while 
CURV  controls  the  number  of  additional  points  to  be  added  in  regions  with  high 
curvature  in  the  concentration  profiles.  Lower  numbers  for  GRAD  and  CURV  mean 
more  mesh  points  will  be  inserted  each  time  the  adaptive  mesh  routine  is  called.  The 
default  values  for  GRAD  and  CURV  are  0.1  and  0.5,  respectively.  However,  it  has  been 
found  that  in  some  cases  a  GRAD  value  of  0.1  could  be  too  restrictive  for  achieving  a 
fully  converged  solution.  This  is  especially  true  near  a  flame  front  where  sharp  changes 
in  radicals'  concentration  create  large  gradients  and  curvatures  in  concentration 
profiles  of  radical  species.  Furthermore,  inserting  more  mesh  points  does  not 
necessarily  improve  the  accuracy  of  the  solution,  but  it  definitely  demands  more 
computing  resources.  Since  Newton  algorithm  is  more  likely  to  converge  on  a  coarse 
mesh,  the  strategy  used  was  to  start  the  computations  with  a  coarse  grid  and  relaxing 
the  GRAD  parameter,  to  say  0.8.  Once  a  solution  has  been  achieved  the  value  of  GRAD 
is  then  reduced  and  the  computations  are  repeated  restarting  from  the  previously 
converged  solution.  This  cycle  is  continued  until  the  required  accuracy  is  achieved. 
Typically,  GRAD  and  CURV  values  of  0.3  and  0.5,  respectively,  have  been  used. 

Convergence  of  the  solution  is  achieved  when  the  residuals  of  the  discretised 
conservation  equations  for  mass  and  energy  are  reduced  to  a  pre-selected  tolerance  and 
when  the  mesh  is  sufficiently  refined  to  meet  specified  gradient  and  curvature  criteria 
of  the  solution  profiles.  Details  on  the  numerical  schemes,  the  adaptive  mesh 
algorithm  and  convergence  criteria  can  be  found  in  Refs[12-15]. 


4.  Model  Parameters 


A  schematic  block  diagram  of  the  structure  of  the  code  is  given  in  Fig.  5.  The  numerical 
model  requires  a  number  of  input  parameters  that  include  gas  phase  kinetics  and 
optionally  surface  reactions,  thermodynamic  and  transport  properties  of  all  species  in 
the  reactions  mechanism.  Also  required  are  the  physical  and  optical  properties  of  the 
porous  medium. 


4.1  Chemical  Kinetics  Model 

All  simulations  have  been  carried  out  for  a  stoichiometric  propane/ air  mixture.  The 
chemistry  of  propane  (C3H8)  has  been  represented  by  a  kinetics  mechanism  with  40 
reactive  species  and  consists  of  93  elementary  reactions.  The  kinetic  scheme  included 
also  NOx  formation  and  rebuming  kinetics  that  have  been  adopted  from  the  GRI  2.1 
mechanism[20].  The  details  of  the  propane  mechanism  are  given  in  Appendix  A. 
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4.2  Numerical  Parameters  and  Physical  Properties 

The  model-input  parameters  include  burner  geometry,  firing  rate,  and  reactants'  mole 
fractions.  The  porous  medium  is  represented  in  the  model  as  a  single  layer,  3mm  in 
thickness  made  of  INCONEL601  fibre  matrix  of  0.4mm  wires.  INCONEL601  is  a  nickel- 
chromium  alloy  ( typical  composition  is  80%  Ni,  16%  Cr,  and  4%  Fe)  which  is  used  in 
many  engineering  applications.  It  has  excellent  mechanical  properties  and  is  resistant 
towards  oxidation  in  air  for  up  1093  °C.  The  total  hemispheric  emissivity  s,  of 
INCONEL601  varies  modestly  with  temperature.  For  example  at  a  temperature  of 
315°C,  the  emissivity  s  is  0.69,  but  it  increases  to  0.82  at  980°C.  Since  the  operating 
temperature  of  the  PRB  is  expected  to  be  in  a  temperature  range  around  900°C,  a 
representative  emissivity  value  of  0.8  has  been  selected  and  used  in  all  calculations. 

For  simplicity  the  computational  domain  used  in  this  study  starts  at  the  upstream  edge 
of  the  porous  medium  ,  ie.  Xi„=xo  =0,  and  extends  to  xout=  10  mm  in  the  downstream 
direction,  see  Fig.  3.  The  physical  and  optical  properties  of  the  porous  medium  such  as 
thermal  conductivity,  emissivity  and  so  forth,  are  also  required.  Table  1  summarises 
the  input  parameters  that  have  been  used  and  kept  unchanged  in  all  the  calculations. 
The  values  of  these  parameters  are  selected  to  closely  represent  those  of  the  experiment 
so  a  consistent  comparison  can  be  made. 


Table  1:  List  of  parameters  that  have  been  kept  unchanged  throughout  all  the  calculations. 


Property  (units) 

Value 

Propane  mole  fraction  (-) 

0.042 

Equivalence  ratio  (-) 

1.0 

Firing  rate  (kW /  m2) 

417. 

Mixture's  inlet  temperature  (K) 

350 

Thermal  conductivity  of  mesh  (W/m.K) 

27.5 

Wire-mesh  thickness  (mm) 

0.4 

Porous  medium  thickness  (mm) 

3.0 

Computational  domain  length  (mm) 

10 

4.3  Sensitivity  Analysis 

To  examine  the  effect  of  variations  (due  to  uncertainty)  in  some  input  values  on  the 
solution,  a  systematic  sensitivity  analysis  has  been  carried  out.  The  gas  and  solid 
temperatures,  concentration  profiles  of  CO,  CO2,  OH,  and  the  net  radiation  flux  have 
been  selected  as  indicators  to  determine  the  solution's  sensitivity  (or  lack  of  it)  to 
variations  in  any  specific  input  parameter.  Sensitivity  analysis  has  been  performed 
against  a  total  of  nine  input  parameters:  three  parameters  that  describe  the  porous 
medium  property  and  structure,  namely  voidage,  tortuosity,  and  area  density,  and  six 
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parameters  that  describe  the  optical  properties  of  the  medium.  These  are  the 
emissivity,  reflectivity,  single  scattering  albedo,  extinction  coefficient,  and  the  forward- 
and  backward-scattering  fractions.  A  list  of  the  various  input  configurations  that  have 
been  examined  is  shown  in  Table  2.  The  approach  used  for  selecting  the  optical 
parameters  shown  in  Table  2,  is  first  to  examine  the  solution's  sensitivity  at  the  extreme 
limits  of  a  certain  parameter,  e.g.  for  the  single  scattering  albedo  (co0)  such  limits  are  0.1 
and  0.9.  Then  the  sensitivity  analysis  is  repeated  for  a  value  that  closely  reflects  the 
actual  property  of  the  porous  material  that  is  used  in  the  experiment.  For  example  a 
single  scattering  albedo  of  0.7  is  typical  for  many  Nickel  based  materials. 


Figure  5:  A  block  diagram  of  the  structure  of  the  ChemRad  -porous  burner  model. 
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Table  2:  Summary  of  the  various  configurations  that  have  been  examined. 

Key:  <f)  =  Void,  y=  Tortuosity,  Aden~  Area  Density (cm2/cm3),  s  -  Emissivity,  go  -  Reflectivity, 
ojo=  Single  Scattering  Albedo,  <je=  Extinction  coefficient  (cm-1),  f&  b  = forward-  and  backward- 
radiation  scattering  fractions,  respectively. 


Configuration 

4> 

Y 

Aden 

e 

Xo 

(Bo 

ae 

/ 

b 

Cl 

0.9 

1.12 

7 

0.8 

0.6 

0.7 

50 

0.5 

0.5 

C2 

0.3 

1.12 

7 

0.8 

0.6 

0.7 

50 

0.5 

0.5 

C3 

0.9 

1.12 

70 

0.8 

0.6 

0.7 

50 

0.5 

0.5 

C4 

0.95 

1.12 

44 

0.8 

0.6 

0.7 

50 

0.5 

0.5 

C5 

0.95 

1.12 

44 

0.8 

0.6 

0.7 

5 

0.5 

0.5 

C6 

0.95 

3.36 

44 

0.8 

0.6 

0.7 

5 

0.5 

0.5 

C7 

0.95 

3.36 

44 

0.8 

0.6 

0.7 

50 

0.5 

0.5 

C8 

0.95 

3.36 

44 

0.8 

0.1 

0.7 

50 

0.5 

0.5 

C9 

0.95 

3.36 

44 

0.8 

0.9 

0.7 

50 

0.5 

0.5 

CIO 

0.95 

3.36 

44 

0.8 

0.9 

0.1 

50 

0.5 

0.5 

Cll 

0.95 

3.36 

44 

0.8 

0.9 

0.9 

50 

0.5 

0.5 

Cl  2 

0.95 

3.36 

44 

0.8 

0.9 

0.9 

50 

0.1 

0.9 

C13 

0.95 

3.36 

44 

0.8 

0.9 

0.9 

50 

0.8 

0.2 

5.  Results  and  Discussion 


Figures  6-17  show  the  modelling  results  that  correspond  to  the  configurations  listed  in 
Table  2.  For  cases  where  the  solution  does  not  change  beyond  a  certain  point,  the 
results  are  displayed  only  over  that  relevant  portion  of  the  computational  domain.  By 
simultaneously  inspecting  the  parameters  shown  in  Table  2,  and  the  corresponding 
figures  it  is  possible  to  identify  the  effect  of  variations  in  one  or  more  of  the  parameters 
on  the  solution.  For  example,  comparing  configurations  Cl  and  C5  allows  the 
combined  effect  of  variation  in  area  density  and  extinction  coefficient  to  be  examined. 

By  comparing  the  results  of  cases  Cl  and  C2,  the  effect  of  void  can  be  identified.  It  is 
shown  that  reducing  the  voidage  of  the  porous  medium  by  a  factor  of  3  causes  an 
increase  in  the  gas  and  solid  temperatures  by  10%,  and  6%,  respectively.  The 
concentration  of  CO2  at  the  exit  plane  of  the  burner  (x=0.3  mm)  has  increased  by 
around  2%.  However,  CO  has  increased  by  25%,  and  OH  by  a  factor  of  2.4.  This  is 
mainly  due  to  the  sensitivity  and  strong  dependence  of  CO  and  OH  concentrations  on 
the  flame  temperature.  The  net  radiation  is  also  increased  (by  a  factor  of  2)  as  expected 
because  reducing  the  voidage  leads  to  an  increase  in  the  residence  time  of  the 
combustion  gases  and  the  contact  surface  of  solid  and  gas,  hence  enhancing  heat 
transfer  within  the  porous  layer. 
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Comparing  Cl  and  C3  shows  the  effect  of  an  increase  in  area  density.  In  this  case  an 
increase  in  area  density  by  a  factor  of  ten  resulted  in  a  reduction  of  the  gas  temperature 
by  approximately  180K,  and  increased  the  solid  temperature  by  approximately  270K. 
However,  while  concentrations  of  CO  and  CO2  have  slightly  increased  and  the  OH 
radicals  have  slightly  reduced,  the  net  radiation  flux  has  increased  by  a  factor  of  four. 
This  is  primarily  due  to  the  increase  in  the  surface  area  between  the  porous  medium 
and  the  hot  combustion  products.  This  has  the  effect  of  enhancing  heat  transfer  from 
the  gas  stream  to  the  solid  hence  lowering  the  gas  temperature  and  increasing  the  solid 
temperature.  A  high  temperature  of  the  solid  combined  with  high  emissivity  causes 
an  enhanced  thermal  radiation  transfer.  The  changes  in  OH  and  CO2  concentrations 
can  probably  be  attributed  to  the  low  gas  temperature,  but  the  cause  for  the  increase  in 
CO  is  not  fully  understood  yet. 


Figure  6:  Axial  distribution  of  gas  temperature  for  configurations  C1-C7  as  listed  in  Table  2. 
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Cl  - (2  - C3  - 04  —A — C5  -*-06  — e— C7 


Figure  7:  Axial  distribution  of  solid  temperature  xmthin  the  porous  medium  for 

configurations  C1-C7  as  listed  in  Table  2. 


Figure  8:  Mole  fraction  profiles  of  CO  for  configurations  C1-C7  as  listed  in  Table  2. 
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Figure  9:  Mole  fraction  profiles  ofCChfor  configurations  C1-C7  as  listed  in  Table  2. 


Figure  10:  Mole  fraction  profiles  of  OH  radicals  for  configurations  C1-C7  as  listed  in 

Table  2. 
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- Cl  - C2  - C3  - CA  —A — C5  -*-05  -&-CI 

Figure  11:  Net  radiation  for  configurations  C1-C7  as  listed  in  Table  2. 


Configurations  C4  and  C5  examine  the  effect  of  reducing  the  extinction  coefficient  by  a 
factor  of  10,  which  interestingly  did  not  record  significant  effect  on  the  results.  In 
general,  however,  changes  in  the  extinction  coefficient  usually  have  profound  effect  on 
the  flame  stability  and  combustion  characteristics  in  a  thick  porous  medium  with  small 
pore  dimensions.  This  is  not  the  case  here  therefore  the  effect  of  the  extinction 
coefficient  is  not  evident.  The  only  noticeable  effect  is  a  small  decrease  in  the  net 
radiation  of  approximately  5%,  as  shown  in  Fig.  11. 

The  effect  of  the  porous  medium  tortuosity  is  examined  by  comparing  configurations 
C5  and  C6.  An  increase  in  the  tortuosity  by  a  factor  of  three  did  not  show  any 
significant  changes  in  the  temperature  or  the  concentration  profiles.  Also,  the 
combined  effect  of  changing  both  the  tortuosity  and  the  extinction  coefficient  ( 
comparing  C4  and  C6)  did  not  show  significant  effect  on  the  results. 

The  effect  of  variations  in  reflectivity,  scattering  albedo  and  the  ratio  between  the 
forward  and  backward  scattering  fractions,  are  shown  in  Figure  12-17.  Comparing 
cases  C8  and  C9,  shows  that  reducing  the  reflectivity  by  a  factor  of  nine,  causes  a 
noticeable  increase  in  the  gas  and  solid  temperatures  by  approximately  110K  (~8%) 
and  300K  (35%),  respectively.  The  net  radiation  also  increases  by  a  factor  of  two, 
however  the  concentration  of  CO,  CO2  and  OH  did  not  change  significantly.  These 
results  indicate  that  heat  transfer  mechanism  by  thermal  radiation  dominates  over 
conduction  heat  transfer. 
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12:  Gas  temperature  profiles  for  configurations  C8-C13  as  listed  in  Table  2. 
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Figure  13:  Solid  temperature  profiles  within  the  porous  layer  for  configurations  C8-C13  as 

listed  in  Table  2. 
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Figure  16:  OH  mole  fraction  profiles  for  configurations  C8-1C13  as  listed  in  Tables  2. 
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Figure  17:  Net  radiation  for  configurations  C8-C13,  as  listed  in  Table  2. 

Comparing  configurations  CIO  and  Cll,  shows  that  a  reduction  in  the  scattering 
albedo  by  a  factor  of  nine  does  not  seem  to  a  have  significant  effect  on  temperatures  or 
concentrations  of  major  species.  The  net  radiation,  however,  is  reduced  by 
approximately  10%.  Configurations  C12  and  C13  examine  the  effect  of  changing  the 
ratio  between  the  forward  and  backward  radiation  scattering  fractions.  By  changing 
the  ratio  from  1:9  to  8:2  (and  even  to  1:1  as  in  most  configurations),  the  only  evident 
effect  is  a  slight  reduction  in  the  net  radiation  of  less  than  5% . 

To  summarise,  the  results  of  the  systematic  sensitivity  analysis  presented  above  show 
that  out  of  the  nine  parameters  examined  four  can  be  considered  influential.  These  are, 
in  order  of  significance,  the  area  density,  void,  reflectivity  and  the  ratio  of  the  forward- 
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to-backward  scattering  fractions.  The  degree  to  which  these  parameters  affect  the 
results,  however,  is  also  influenced  by  the  thickness  of  the  porous  layer  and  by  catalytic 
surface  reactions  (which  have  not  been  included  in  this  study).  Therefore,  it  is  essential 
to  recognise  that  this  sensitivity  analysis  is  used  in  the  numerical  sense  to  estimate  the 
effect  of  variations  in  input  parameters  on  the  solution.  This  does  not  necessarily  reflect 
actual  physical  conditions.  That  is,  in  reality  there  is  an  inherent  coupling  between 
various  optical  properties,  which  implies  that  it  is  not  always  possible  to  alter  one 
parameter  while  keeping  the  other  unchanged.  Fu  et  al  [21],  for  example  showed  that 
for  cellular  ceramic  the  single  scattering  albedo  increases  with  the  porosity.  They  also 
showed  that  an  increase  in  the  material  reflectivity  leads  to  an  increase  in  the  scattering 
albedo.  Therefore,  care  should  be  exercised  when  assigning  optical  properties  in  the 
model  to  predict  the  combustion  and  radiative  performance  of  actual  porous  materials. 


6.  Concluding  Remarks 

A  one-dimensional  numerical  model,  ChemRad,  for  predicting  combustion  and  heat 
transfer  characteristics  in  porous  media  has  been  introduced.  The  model  accounts  for 
detailed  chemical  kinetics  and  uses  a  two-flux  radiation  model  in  the  energy  equation. 
The  capabilities  of  the  model  to  incorporate  detailed  gas  and  surface  chemistry 
schemes,  handle  multiple  layers  of  different  porous  materials,  and  couple  the 
chemistry-radiation  provide  a  flexible  tool  that  can  be  used  to  optimise  the 
performance  of  porous  radiant  burners. 

The  main  source  for  potential  inaccuracies  in  the  numerical  solution  arises  from 
uncertainties  in  the  values  of  certain  material  properties,  the  greatest  being  the  optical 
and  the  microstructural  properties  of  the  porous  medium.  A  sensitivity  analysis  has 
been  performed  to  quantify  the  extent  of  the  solution's  sensitivity  to  nine  input 
parameters.  It  has  been  found  that  the  most  influential  input  parameters  in  order  of 
significance,  are  the  area  density,  voidage  and  reflectivity  of  the  porous  material,  and 
the  ratio  between  the  forward-  and  backward  scattering  of  radiation. 

The  robustness  of  the  ChemRad  model  has  been  successfully  demonstrated  in  a 
separate  study [5],  in  which  numerical  predictions  have  been  compared  favourably 
with  experimental  results. 
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Appendix  A:  Chemical  Kinetics  Mechanism 


A  high  temperature  kinetics  mechanism  for  propane  reactions.  The  mechanism 
includes  also  N-O  reactions  from  the  GRI-2.1  mechanism. 
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(K  =  A*Tb 
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